Resonant Inelastic X-ray Scattering (RIXS) Spectra for Ladder Cuprates 
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The ladder compound Sri4Cu2404i is of interest both as a quasi-one-dimensional analog of the 
superconducting cuprates and as a superconductor in its own right when Sr is substituted by Ca. In 
order to model resonant inelastic x-ray scattering (RIXS) spectra for this compound, we investigate 
the simpler SrCu2 03 system in which the crystal structure contains very similar ladder planes. We 
approximate the LDA dispersion of SrCu2 03 by a Cu only two-band tight-binding model. Strong 
correlation effects are incorporated by assuming an anti-ferromagnetic ground state. The available 
angle-resolved photoemission (ARPES) and RIXS data on the ladder compound are found to be in 
reasonable accord with our theoretical predictions. 

PACS numbers: 74.72-h, 75.50.Ee, 78.70.Ck 
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I. INTRODUCTION 



Resonant inelastic x-ray scattering (RIXS) is a second- 
order optical process in which there is a coherent ab- 
sorption and emission of X-rays in resonance with elec- 
tronic excitationsii RIXS can probe charge excitations 
extending to fairly high energies of up to ~ 8 eV. This 
allows the analysis of electronic states over a wide energy 
range, including electron correlation effects originating 
from strong electron-electron Coulomb repulsion, provid- 
ing thus a powerful tool for investigating Mott physics in 
solids. 

The chain-ladder compound Sri4Cu2404i exhibits very 
interesting magnetic, transport and properties. It has 
attracted wide attention due to the discovery of a su- 
perconducting phase in highly Ca-doped samples at high 
pressure^i and charge order of the doped ladder—. The 
compound possesses an incommensurate layered struc- 
ture consisting of alternating layers of sublattices involv- 
ing Cu02 chains and CU2O3 ladders. The supercon- 
ductivity arises on the ladders, making them a quasi- 
one-dimensional analog of the cuprates. Very recently, 
K-edge RIXS data on the ladder compound has been 
reported^, providing motivation for undertaking corre- 
sponding theoretical modeling of the spectra. Here, we 
attempt to do so by considering the simpler analog com- 
pound SrCu203. This should be a good approximation 
since interlayer coupling in Sri4Cu2404i is negligibleS, 
and both compounds have very similar ladder planes 
with similar hopping parameters^. Specifically, we ob- 
tain K-edge RIXS spectra within a mean field approach 
for momentum transfer along as well as perpendicular to 
the direction of the ladders. A two-band Cu-only tight- 
binding model is used in which strong correlation effects 
are incorporated by treating an antiferromagnetic (AFM) 
ground state. 



II. ELECTRONIC STRUCTURE AND THE 
TWO BAND MODEL 



The spin-ladder compound SrCu2 03 possesses the or- 
thorhombic structure with space group Cmmm in which 
CU2O3 planes are stacked with Sr atoms sandwiched be- 
tween these planes. * Fig. 1 shows the detailed arrange- 
ment of Cu and O atoms in the CU2O3 planes. This so- 
called 'trellis structure' involves Cu-0 ladders where suc- 
cessive ladders are seen to be offset by half a unit cell. We 
obtained the band structure of SrCu2 03 self-consistently 
using a full-potential, all electron scheme within the lo- 
cal density approximation (LDA)^ii2,. The first principles 
bands were fitted by a 2-band tight-binding (TB) model 
in the vicinity of the Fermi energy, and provided the ba- 
sis for RIXS computations presented in this study. Fig. 
2 shows the first-principles as well as the TB bands along 
several high symmetry lines in the Brillouin zone (BZ). 
There are seen to be only two bands around the Fermi 
energy, which display large dispersion along the ladder di- 
rection F-Z, and a relatively smaller dispersion along the 
perpendicular F-X direction. In the first-principles band 
structure, both these bands are dominated by states of 
Cu dx^-y^ character whose weight is given by the color 
bar on the right hand side of Fig. 2. 
In making a TB fit to the two aformentioned LDA bands 
near the Fermi energy, we have adapted a Cu-only 2- 
band model suggested in Ref. 7. The detailed form of 
the TB Hamiltonian is discussed in the Appendix. The 
TB bands are seen from Fig. 2 to provide a good fit to the 
LDA bands near the Fermi energy. Our values of various 
parameters, i.e. the on-site energy eo and the hopping 
parameters ti—tg, are seen from Table 1 to be in reason- 
able accord with those of Ref. 0. The meaning of specific 
overlap terms involved in defining ti — tg is clarified by the 
red arrows in Fig. 1 . The present TB model includes not 
only the nearest neighbor hopping terms ti —t^, but also 
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FIG. 1: (color online) Arrangement of Cu and O atoms in 
the CU2O3 planes in the form of a series of ladders offset 
by half a unit cell along the c-axis. Unfilled circles give the 
location of out-of-plane Sr atoms. The green parallelepiped 
marks a primitive unit cell around which the orientation of 
Cu d^2_y2, O Px, and O py orbitals is shown. Red arrows 
give the specific hopping parameters used in two-band model 
fits to the band structure. The AFM ordering of spins in the 
ladders is depicted by the green arrows. 



the longer range hoppings — tg. Interestingly, we find 
that the inter-ladder dispersion (i.e. along Z-A in Fig. 2) 
cannot be fitted well using only a nearest-neighbor hop- 
ping model. Table 1 shows that the intra-ladder hopping 
parameters (ti, t2 and t^) are generally larger than the 
inter-ladder terms such as t^, and tg- This can be un- 
derstood with reference to Fig. 1 where orientation of the 
Cu-(ia:2_y2 and 0-px and Py orbitals is sketched on a few 
sites. An intra-ladder Cu-O-Cu path with a bond angle 
of 180° (e.g. Cu2-0i-Cui) will be expected to provide 
a larger orbital overlap than an inter-ladder path with 
a 90° bond (e.g. Cu4-0i-Cui). The glide symmetry of 
SrCu203 leads to some dispersion anomalies, including 
extra degeneracies at the zone boundaries and an ap- 
parent 47r periodicity of the dispersion along the ladder. 
A 27r symmetry can be effectively restored by including 
different cuts along kx as shown in Fig. 3(a). Similar 
anomalies in c-axis dispersion due to a glide symmetry 
are also found in Bi2Sr2CaCu208 (Bi2212)". The trel- 
lis compound shows the presence of short-range spin or- 
der with a spin gap which is consistent with theoretical 
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FIG. 2: (color online) First-principles LDA band structure of 
SrCu203 along several high symmetry lines in the irreducible 
Brillouin zone is shown by dots of various colors. The colors 
of the dots code the weight of Cu d^2_y2 character in the as- 
sociated wavefunctions as indicated by the color-bar on the 
right hand side of the figure. Solid red lines give the 2-band 
tight-binding model fit to the LDA bands near the Fermi en- 
ergy. 



TABLE I: TB parameters for 2-band model. 



Parameter 



This work 



Ref. 7 



eo 
t-L 

t2 

t3 

t4 
t5 
U 

ts 
tg 



-0.0350(eV) 

0.5650 

0.3800 

0.0400 

0.0520 
-0.1200 

0.0700 

0.0750 

0.0057 
-0.0115 



-0.0450(eV) 

0.5650 

0.3950 

0.0400 

0.0500 
-0.1150 

0.0400 

0.0750 

0.0050 
-0.0200 



predictions^^. Due to the 180° Cui-0i-Cu2 bonds, the 
spins are strongly coupled antiferromagnetically on the 
legs and the rungs of the ladders as indicated by green ar- 
rows in Fig. l.^'^ However, the displacement of successive 
ladders with respect to each other frustrates the develop- 
ment of long range AFM order. One nevertheless expects 
the electronic system to experience significant AFM fluc- 
tuations, which are presumably sufficient to impose an 
underlying dispersion characteristic of the AFM order. 
In this spirit, we have approximated the correlation ef- 
fects within a Hartree-Fock model of an itinerant AFM, 





FIG. 3: (color online) (a) Paramagnetic, and (b) AFM dis- 
persions within the 2-band TB model as a function of fcz at 
kx ~ (red lines) and kx = 2n/a (blue lines). Experimental 
ARPES datapoints are shown as green dots.'i^ 



as in the planar cupratesi^. Taking the on-site energy 
to be U= 3.3 eV ~ 6t, the magnetization was computed 
self-consistently to be to=0.43. The AFM Hamiltonian is 
given in the Appendix and the resulting dispersions are 
shown in Fig. 3(b). Comparison with the paramagnetic 
solution shows that a large gap of ~2.3 eV opens up be- 
tween the upper (UMB) and the lower magnetic bands 
(LMB). The theoretical LMBs display the characteristic 
backfolding near fc^ = 7r/2c, which is in accord with the 
experimentally observed dispersion (green dots in Fig. 
3(b) via ARPES^^, and is reminiscent of a similar effect 
in the insulating planar cuprates. 



III. RIXS SPECTRA 



Our computations of the the K-edge RIXS cross section 
for the Cu Is 4p core level excitation are based on the 
expressioni^iii 



where 



I{q,uj,uj,) = {2t:)^N\W{uj,uj,)\'^ 
xE„.kX;-'(q,k,^)|Mzj(k)|2 



X;,-, (q, k, Lu) = 5{co + E, (k) - Ey (k -f q)) 
xnj(k)[l -njv(q + k)], 



(1) 



(2) 



rij [k) is the electron occupation of the band and Ej (k) 
is the corresponding energy dispersion obtained by self- 
consistently solving the two-band AFM Hamiltonian (see 
Appendix), and 



FIG. 4: (color online) Color plots show computed RIXS spec- 
tra from SrCu2 03 for momentum transfer along (a): The lad- 
der direction and (b): The direction perpendicular to 
the ladders. Theoretical spectra are not resolution broadened 
in order to highlight their intrinsic structure. Yellow vertical 
line marks the momentum transfer value where the detailed 
spectra, including features marked A-C, are considered in Fig. 
5. 



z5(wi.ki)_D(w/,k2) ■ 



(3) 



Here, D(uj, k) = Lo+eis—Sip{ki)+iTis, 7 is the matrix el- 
ement for scattering from Is to 4p, and Vd is the core- hole 
potential in 3c? level. LUi^LUf) and qi{qf) denote the initial 
(final) energy and momentum, respectively, of the pho- 
ton, and oj = LtJi — ijj f and q — Qi — qf give the energy and 
momentum transferred in the scattering process. Since 
Cu Is is a core state, the associated energy band eis{k) 
is assumed dispersionless. The Cu 4p band dispersion 
e4p{k) is modeled by a 2D-TB model with nearest neigh- 
bor hopping. Tis is the decay rate of core hole taken to 
be 0.8 eV. The matrix element Mij associated with the 
interaction between the core hole and 3d levels around 
the Fermi energy is 



la 



(k 



(4) 

in terms of the eigenvectors X^^ of the AFM Hamil- 
tonian, where a denotes electron spin and I an orbital 
index, ai = Vi/Vd, where Vi is the Coulomb interaction 
between a core hole and an electron on atom I separated 
by a distance i?/. Here we approximate the vertex 
correction A — > 5a^a'. 

Fig. 4 shows RIXS spectra computed within the 2- 
band AFM model in the form of a color plot for mo- 
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FIG. 5: (color online) (a) Lowest and uppermost tight-binding 
AFM bands in SrCu203 at = 0. (b) Red line gives partial 
contribution to the RIXS spectrum at = 7r/2c (marked 
by yellow vertical line in Fig. 4(a)) from transitions between 
the two bands in (a). Blue line gives the corresponding total 
RIXS spectrum. Features A-C are discussed in the text. 



mentum transfer along as well as perpendicular to the 
direction of the ladders. The d-band spectra have not 
been broadened in order to emphasize the presence of 
considerable intrinsic structure in the spectra, despite 
the large broadening F associated with the short core 
hole lifetime. Insight into the nature of these spectra 
can be obtained by examining expressions 1-4 on which 
the computations are based. The enhancement factor 
WiiOi^LOf) is found to vary relatively slowly with en- 
ergy due to the large 4p bandwidth and the substan- 
tial damping of the core hole given by Fig. There- 
fore, spectral shapes are controlled effectively by the 
term -f E^i^) - Eyi:^ + q)) x nj(k)[l - nj/(q + 

k)] I X/^([j,k)X;{,(cj,k+q) |2. Note that this in- 
volves not only the joint density of states (JDOS) factor, 
nj(k)[l -nj/(q-hk)]5(w-t-£;j(k) -i;j/(k-|-q)), but also 
the partial electron occupancy of the filled band given 
by I Xl^[uj^\f.) p and the partial electron occupancy 

I Xj-'^, (w, k+q) p of the empty band. A large contri- 
bution thus results when JDOS connects band extrema, 
leading to resonant peaks in the RIXS cross-section. 
Fig. 5 considers the spectrum at — it /2c in greater 
detail (i.e. corresponding to the vertical yellow line in 
Fig. 4(a)). The blue curve in Fig. 5(b) gives the total 
RIXS cross-section, which of course involves contribu- 
tions from all allowed transitions from either of the two 
unfilled bands to one of the two empty bands in the AFM 
band structure of Fig. 3(b). The red curve in Fig. 5(b) 
gives the partial contribution to the spectrum from just 
the pair of bands shown in Fig. 5(a), i.e. the lowest oc- 
cupied and the highest unoccupied band. In particular. 
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FIG. 6: (color online) Comparison of unbroadened (blue 
curves) and broadened (black curves) RIXS spectra with the 
corresponding experimental data (red dots) taken from Ref.Q. 
Left hand side panels are for gyc at various values, while 
the right hand side panels are for q^c at various q^ values as 
indicated. Spectra are normalized as shown. 



peak C around 3.5 eV in Fig. 5(b) arises from transi- 
tions in (a) marked by the green arrow, while peak B 
has its origin in the transitions given by the orange ar- 
row. [Note that the horizontal shift in the direction of 
the arrows in Fig. 5(a) is the momentum transfer vector, 
while the vertical displacement is the energy transferred 
in the scattering process.] Other spectral details can be 
analyzed in a similar manner and associated with specific 
transitions by examining the partial contributions from 
various pairs of bands. In particular, the resonant peak 
A around 2.5 eV in Fig. 5(b) results from transitions 
between the uppermost filled band and the lowest empty 
band in Fig. 3(a). Along these lines, the intense feature 
around 3 eV in the q± spectra of Fig. 4(b) is found to be 
associated with transitions between the uppermost filled 
band and the lowest empty band (as a function of g^). 

Fig. 6 compares our theoretical spectra with the avail- 
able experimental RIXS data of Ref. 12 on the ladder 
compound. Left hand side panels are for momentum 
transfer along the ladder direction (i.e. gyc) with vary- 
ing from to tt/c, while the right hand side panels are 
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for q^c with varying over the range O-n/a. The un- 
broadened theoretical spectra (bhic hnes) require a sub- 
stantial broadening for a meaningful comparison with the 
data (red dots). Accordingly, we have applied a com- 
bined Gaussian and Lorentzian broadening to the com- 
puted spectra to obtain the broadened theoretical spec- 
tra in Fig. 6 (black lines). The Gaussian broadening 
is taken as the nominal experimental resolution of 120 
meV— . The residual broadening, which reflects presum- 
ably lifetime effects not accounted for in our computa- 
tions, is modeled via a Lorentzian with half- width- at- 
half- maximum of rd=l.l eV for q \\ c and of rd=1.4 eV 
for q _L (Mi^. Although Figs. 4-6 show that the RIXS 
spectra intrisically contain considerable information con- 
cerning the charge excitations and their momentum de- 
pendencies in the ladder compound, much of this struc- 
ture is seen to be lost in the broadened spectra. Never- 
theless, for both II c as well as g ± c, the experimental 
spectra are in reasonable accord with the broadened the- 
ory, some discrepancies in shape and fine structure in 
the experimental data notwithstanding. In particular, 
the spectra show a dispersion of ^ O.SeV along the lad- 
ders (qz) and negligible dispersion perpendicular to the 
ladders (q-x)- While we have compared our calculations 
to the data of Wray et al.^, the undoped data of Ishii et 
al^ show very similar broadening and dispersion, with a 
slightly larger gap. 

IV. CONCLUSIONS 

We have presented a two-band model and the asso- 
ciated K-edge RIXS spectra for the ladder compound 
SrCu203 as a way of capturing the physics of the more 
complex ladder compound Sri4Cu2404i. RIXS spectra 
are considered for momentum transfer along as well as 
perpendicular to the direction of the ladders. Our anal- 
ysis indicates that the available ARPES and RIXS data 



are consistent within experimental resolution with the 
presence of strong antiferromagnetic correlations in the 
system, which are modeled in our study by introducing 
an ordered AFM state. Notably, our calculations do not 
require any significant renormalization of LDA-based 
band paramters in the ladder compound, and moreover, 
our value of the effective Hubbard U is similar to that 
found in the insulating planar cuprates^^. 
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Appendix: AFM and Paramagnetic Hamiltonians 
and Dispersions 

The AFM Hamiltonian for the two band model can be 
written as 

H = £od+dj + E<.,,> {d+d, + d,d+)+ 

where ridj — dj'dj, Eq is the on site energy, t^ 's are the 
hopping parameters (see Fig. 1), and U is Hubbard U. 
The AFM ordering of spins is shown in Fig. 1 by green 
arrows. The Hartree-Fock decomposition of the Hubbard 
term in the Hamiltonian is given by 

rid^f < Udji > + < nd,T > "rfji- 

Urn = y(< "-djT > ^ < "rfii >) = ^"^ 

< n > = < Ud^i > + < Ud^i > 
We can now rewrite the AFM Hamiltonian as 



where 



Hn = H44 = A + U,n, W22 = W33 = A - C/„, (7) 

Ti-u — — 12 — exp(ifcj,/2) — exp(3iA;2/2) — 2^8 cos(2fcz), 

Wis = -2ticos{kz) -2Ucos{{kz - kx)/2) -2t9Cosii3k, ~ k^)/2), 

T-iii = -2t5Cos(fc2) - t3exp(-i(fc^ - fcj.)/2) - i6exp(-i(3A;^ + fc^)/2), 

H23 = -2t5Cos(fc2) --t3exp(i(fc2 - fc^)/2) - t6exp(i(3/c2 - fc^)/2), 

H24 = -tiCOs(fc^)-2t4Cos((fc^ + fc:,)/2)-2t9Cos((3fc^4-fc^)/2)), 



A = £0 - 2t7Cos(2fc^) - 2t4Cos((fc2 + fca;)/2) - 2t9COs((3fc^ - fc^)/2)). 

For the the paramagnetic case Um is equal to zero and the 4x4 AFM Hamiltonian is reduced to a 2 x 2 form with 
matrix elements 

Hii = H22 = eo - 2ii cos(fc^) - 2t7 cos(2fc2) - At^ cos{kx/2) cos(A:^/2) - Atg cos{kx/2) cos{3kz/2), (8) 
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Hi2 = — <2 — 2t5 cos{kz) — 2ts cos(2fcz) — 2t3 exp{ikx/2) cos{kz/2) ~ 2te exp{ikx/2) cos{3kz/2), 
The resulting dispersion ia^ 

e±(k) = eo + cos(fc,/2) + (ei,i(fc^) cos(fc + x/2) ± {e^Akz? + e±,4(fcz)' + 2ei,3(/c,)ei,3(fc.) cos(k,/2)fl^ (9) 

where 

e||(/c^) = -2tiCos(fc2) - 2t7Cos(2/c^), (10) 
ej-,i(^^) = -4t4Cos(fc2/2) - 4t9COs(3fc^/2), 
ej-,3(^^) ^2 + 2i5COs(fc2) + 2i8COs(2/c^), 
C-A^Ak.) = 2i3Cos(fcj2) + 2t6COs(3A:j2), 
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